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1. Introduction 
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Abstract. Motivated by sample path decomposition of the stationary continuous state branch- 
ing process with immigration, a general population model is considered using the idea of im- 
mortal individual. We compute the joint distribution of the random variables: the time to the 
most recent common ancestor (MRCA), the size of the current population and the size of the 
population just before MRCA. We obtain the bottleneck effect as well. The distribution of the 
Z^ , number of the oldest families is also established. The results generalize those in the recent paper 

by Chen and Delmas [8]. 

< 

•^ ■ Continuous state branching processes (CB-processes) are non-negative real-valued Markov 

1 -Q , processes first introduced by Jirina [T7] to model the evolution of large populations of small 

C^ I particles. Continuous state branching processes with immigration (CBI-processes) are gener- 

alizations of those describing the situation where immigrants may come from outer sources, 
see e.g. Kawazu and Watanabe [18]. It is shown in Lamperti [22] that a CB-process can be 
obtained as the scaling limit of a sequence of Galton- Watson processes; see also [5] [H [23]. A 
genealogical tree is naturally associated with the Galton- Watson process. This has given birth 
, to the continuum random tree theory first introduced by Aldous [H 1^ to code the genealogy of 

Q ■ the CB-process. Duquesne and Le Gall [S] further developed the continuum Levy tree to give 

^SJ , the complete description of the genealogy of the CB-process in (sub)-critical case. Kingman 

,— b ' has initiated the study of the coalescent process in 1982 in his famous papers [W\ I2UJ . Then 

^^ ■ coalescents with multiple collisions, also known as A-coalescents, were first introduced and stud- 

^^ I led independently by Pitman [26] and by Sagitov [27]. Recently some authors have studied the 

coalescent process associated with branching processes, see e.g. Lambert [21] on coalescent time, 
Evans and Ralph [11] on the dynamics of the time to the most recent common ancestor (MRCA) , 
^ I Chen and Delmas [8] on MRCA on some special stationary CBI- process, and Berestycki et al.[7] 

H ' on the coupling between A-coalescents and branching processes. 

C^ ■ 

This paper is motivated by Chen and Delmas [8]. The model considered here is a direct 

extension of [8]. We will use some notations and definitions in that paper and consider the general 

CBI-process here instead. The fact that the CBI-process may have a non-trivial stationary 

distribution makes it a more interesting model to be considered here than the CB-process since 

for the CB-process either the population becomes extinct or blows up with positive probability. 

We consider a (sub)-critical CBI-process Y = {Yt,t > 0) with branching mechanism ip given by 

(j2.1|) and immigration mechanism F given by (j2.7p . Our main interest is in presenting a further 

model of random size varying population and exhibiting some interesting properties. Afterwards 

we will give some properties of the coalescent tree. 
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We consider the stationary CBI-process defined on the real hne Z = {Zt,t € M). In order for 
the time to MRCA to be finite, we assume condition (Al): 



~ dz 
1 i'iz) 



In order for Zt to be finite, we shall assume condition (A2): 

"^ F{z) 



^{z] 



■ dz < oo, for some A > 0. 



Using the look-down construction for the population with constant sizes, we represent the 
process Z by means of the picture of an immortal individual which gives birth to independent 
populations. We first give some notations. For fixed time t = (indeed we can choose any time 
by stationarity) , we denote by A the time to the MRCA of the population living at time 0, 
Z = Z/_y^\_ the size of the population just before MRCA, and Z^ the size of the population at 
time which has been generated by the immortal individual over the time interval {—A, 0) and 
Z'^ = Zq — Z^ the size of the population at time generated by the immortal individual at time 
—A. We will see that conditionally on A, the random variables Z , Z and Z are independent, 
and the joint distribution of the random variables is also considered. We also obtain the result 
that the size of the population just before MRCA is stochastically smaller than that of the 
population at the current time, that is the bottleneck effect. 

Let A^ + 1 represent the number of individuals involved in the last coalescent event of the 
genealogical tree. We present the joint distribution of A, N and Zq. Using the measured rooted 
real tree formulation of the genealogy of the stationary CBI-process developed in (2j, we give 
the asymptotic for the number of ancestors. 

We will give the transition probabilities of the MRCA age process {At ,t S M) , which has 
been studied by Evans and Ralph in |llj for the CB-process conditioned on non-extinction. We 
generalize it to the general case with the similar lines as their proof. In the end we study the 
zero set of the CBI-process as well, which is a stationary regenerative set. Foucart and Bravo 
|12j have studied the CBI case on the positive half line. The stationary case is a bit different 
as the subordinator is not naturally associated with the regenerative set; see j28j and [16] for 
details. For this situation see also 1151. 



This paper is organized as follows. We first recall some well-known results on the CB-process 
and CBI-process in Section 2. The family and clan decomposition of the CBI-process are then 
introduced in Section 3. We will give the condition for the existence of the stationary CBI- 
process, determine the joint distribution of A, Z , Z^ , Z'^ and prove the bottleneck effect in 
Section 4, that is Z is stochastically smaller than Zq. In Section 5 the distribution of the 
number of individuals involved in the last coalescent event A^ is computed. In the latter part 
of Section 6 we will introduce the genealogy of CB-process using continuum random Levy trees. 
Then the asymptotic for the number of ancestors is given. In Section 7 we give the transition 
probabilities of the MRCA age process and the properties of the zero set. 



2. CB-PROCESS AND CBI-PROCESS 

We recall some well-known results on CB-process and CBI-process derived from Li \24:\ [25] . 
We consider a (sub)-critical branching mechanism ^: 

/•oo 

(2.1) '4){z) = hz + cz'^ + j (e-^"-l-Fzn)m(dii), z > 0, 

JO 
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where b = ^ (0+) > 0, c > are constants and (n A u'^) m{du) is a finite measure on (0, oo). We 
will consider the non-trivial case, that is, assumption (A3): 



either c > or / um{du) = oo. 
^(0,1) 



(0,1) 

There exists a M+-valued strong Markov process X = {Xt,t > 0) called continuous state branch- 
ing process (CB-process) with branching mechanism ip whose distribution is characterized by 
its Laplace transform 

(2.2) E^.[e-^-'^*] =e-^'^*(^\ 

where E^^ means that Xq = x and the function ft (A) is the unique non-negative solution of the 
backward equation 

(2.3) ; §iM>^) = -HM>^)), t > 0, A > 0, 

vo{X) = A, A > 0, 

The CB-process has a canonical Feller realization. Let P^; be the law of such a CB-process 
started at mass x > 0. Moreover, X has no fixed discontinuities. The probability measure 
Qt{x, •) is infinitely divisible and under condition (A3), vt{X) can be expressed canonically as 

/■oo 

vt{X)= (l-e-^")/i(dn), t>0,A>0, 

Jo 

where ult{du) is a finite measure on (0, cxo); see Theorem 3.10 in p4j. The Markov property of 
X implies that for any X,s,t > 0, 

(2.4) vt+s{X) = vt{vsiX)). 
We also have the forward differential equation 

(2.5) ^^vt{X) = -^P{X)^vt{X). 

Let C = inf{s > 0,Xs = 0} be the extinction time of X and c{t) = limA-i.oo vt{X). Under (Al), 
c{t) > is finite. We have by ([23]) that 

(2.6) Vs{c{t)) = c{t + s). 
We consider an immigration mechanism F: 

/>oo 

(2.7) F{z) = l3z+ (l-e-^")n(du), z>0, 



where /3 > is a constant and (1 A u) n[du) is a finite measure on (0, oo). Then there exists a 
strong Markov process Y = {Yt,t > 0) called continuous-state branching process with immigra- 
tion (CBLprocess) with branching mechanism ip and immigration mechanism F defined on M-|- 
with Laplace transform given by 

(2.8) ^x[e~^^'] =e-^'^*W-i'o^("-W)'^^ 

where E^, means that Yq = x. We also denote P the corresponding probability measure. 
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3. Sample path decomposition 

In this section we will recall some results from Li |241 [25] . We will give the clan and family 
decomposition of the CBI-process. 

The process X is infinitely divisible. It is well known that there exists a canonical measure 
(we also call it excursion law) Qo on the space D of Cadlag functions on [0, oo) with Skorokhod 
topology. Notice that Qo({X,Xo+ / 0}) = 0. 

We can give a reconstruction of the sample paths of the CB-process by means of the excursion 
law. Let X >0 and let N{dX) = Xlie/ ^x^ (dX) be a Poisson random measure on D with intensity 
xQQ{dX). We define the process {X^,t > 0) by 



(3.1) 



Xo 



X[= I XtN{dX), t>0. 
Jd 



Then X is a realization of the CB-process X. We will not distinguish X from X. For the 
proof see Li |241 Theorem 8.24] or |25[ Theorem 2.4.2]. As one can see ()3.ip is equivalent to the 
well-known decomposition as follows: If N{dx,dX) = Ylii=j St^.x'^){dx,dX) is a Poisson point 
measure on R+ x D with intensity l[Q^(y2-j{x)dxQo{dX), then ^^^^ 1|3,.<^|X* is distributed as 
X under P^,.. Further we have for A > 0, 

Qo(l - e-^^') = lim -E,[l - e"^^'] = vt{X), 

x^O x 

and 

c{t) = Qo(C > t) = Qo(^t > 0). 
We will put Xt = for t < 0. 

Now we will introduce the family decomposition of the CBI-process. We consider the following 
Poisson point measures. 

(1) Let NQ{dr,dt) = "^i^j S(^rt,ti){dr,dt) be a Poisson point measure on (0,oo) x R with 
intensity n{dr)dt. 

(2) Conditionally on A'o, let {Ni^i,i € /) be independent Poisson point measures with inten- 
sity ri6ti{dt)Qo{dX) , where Ni^i{dt, dX) = J^jej^ i ^{t-,x^){dt, dX). Note that for all j G 
Ji,i, we have tj = U. We set Ji = UiG/^i.«' ^"^ Ni{dt,dX) = J2jeJi ^{tj,x^)idt,dX). 

(3) Let N2{dt,dX) = '^j^j 6u.xi){dt,dX) be a Poisson point measure with intensity f3dt 
Qo(dX) independent of Nq,Ni. 

We set J' = Ji\]J2- We shall call X^ a family and tj its birth place for j G J. We will 
consider the process {Y^ ,t > 0) and its stationary version {Zt,t G M). They are usually called 
family decomposition of the CBI-process defined as follows: 

(3.2) Yl= Y. ^t-t,, Z, = Y,Xl_,. 

jeJ,tj>o jej 

Putting this another way we can deduce that it corresponds to a special immigration process 
shown in Corollary 3.4.2 in |25j . 

For i G / denote X* = X^jgj--^ . X^ and I = I\J J2- The random measure 

(3.3) N^idt.dX) = Y,Siu,x^){dt,dX) 
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is a Poisson point measure with intensity dtfi{dX), where fi is given by 
(3.4) n{dX) = /3QoidX) + [ n{dx)P^{dX). 

J(0,oo) 

It corresponds to the entrance law {Ht,t > 0) in Li [25] given by 

/>oo 

Ht = (3lt+ / n{dx)Qt{x,-). 
Jo 

We shall call X^ with i G X a clan and ti its birth place. For j G J^2, X^ is a clan and a family. 
Then we have 



(3.5) 



iex,t,>o iex 



Y with this representation is just the sample path decomposition of the CBI-process. We shall 
call this the clan decomposition of y. Y is a version of Y and Z is the stationary version of Y. 
Usually the family decomposition is more precise than the clan decomposition. 

We give an interpretation of Z in population terms. At time t, Zt corresponds to the size of 
the population generated by an immortal individual giving birth at rate /3 with sizes evolving 
independently as X under Qq and at rate 1 with intensity n{dx) with initial size x which evolve 
independently as X under P^;. We first give a lemma on the family representation. 

Lemma 3.1. Let f be a non-negative measurable function. We have 



(3.6) 



E 



'E,6j/fe,^^) 



7.F(Qo(l-c-/('.^)))dt 



Proof. Due to the independence of the Poisson random measures and the exponential formula, 
we have 

E 



e-E,ej/fe.^^)" 


= E 


'^~j:,ej,m<x^)' 


E 


e 


E.e^^/fe'^')" 






= E 


■g-E,e,E,ej,,,m,^^r 


g-/3 4Qo(l-e-/(*.^))dt 




= E 


e-E.einQo(l-e-^(*-^)) 


" g-/3/^Qo(l-c-/('>^))dt 



^-J^FiQoil-e-f(^-^)))dt 



D 



The necessary and sufficient condition for which the CBI-process has a stationary version is 
that (A2) holds, see Theorem 3.20 in |24j . If (A2) holds, then Xt converges in distribution to 
Xoo as t — )• oo, with the distribution of X^o characterized by its Laplace transform 



-AXo 



(3.7) E[e" 

Then with (A2) in force, Z defined by (|3.2p and (j3.5p is the stationary version of Y. 

Corollary 3.1. Assume that (A2) holds. We have for A > 0, t G M, 

(3.8) 

In particular, we have 

(3.9) E[Zt] 



E[Ziexp(-AZ,)] = :^E[e"^^* 
'0(A) 

F'(0) 



V^'(o)- 
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Proof. We can see from ()3.7p that 

/•oo 

E[Zt expi- XZt)] =E[e-^^']dx / F(t;,(A))ds. 

JO 

Using the forward equation ()2.5p we can deduce that 

/•oo /'OO 

Jo Jo 

1 /-"^ „/ „ ... , F(A) 



The second part is obvious. 

In the fohowing we wih always suppose that (Al), (A2) and (A3) are in force. 



D 



4. Time to MRCA and the population sizes 

With the decomposition procedure in force, we wih foHow the steps of Chen and Delmas [8]. 
We consider the coalescence of the genealogy at a fixed time to. We may as well assume that 
to = because of stationarity. There are infinitely many number of clans contributing to the 
population at time 0. We can further prove that there are only finite number of clans born 
before time a and still alive at time 0. Only one oldest clan is expected to be still alive at time 
0. 

First we will give the notations using the decomposition. —A is the birth time of the unique 
oldest clan at time {A is also the time to the most recent common ancestor (TMRCA) of the 
population at time 0) given by ^ := — infjtj < 0, ^it, > 0, i E X}; Z^ is the population size of 
this clan at time 0, i.e. Z*^ := ^i^,, if j4 = — tj; The size of all the clans alive at time with 
birth time in (—^4,0) is given by Z^ := Zq — Z'^, and the size of the population just before the 



MRCA is given by Z^ := Z, 



i~Ay 



E^exXt 



A-u'^lUK-A}- 



Theorem 4.1. Let f -.R^l 



be a measurable function. For A,7,r? > 0, we have 



dtf{t){Fic{t))-F{vtir^)))exp 



t l-OO 

F{vs{-i))ds- / F{vs{\ + c{t)))ds]. 
Jo 



Proof. We have 












^^^-XZ^-,Z'-,ZO ^(^^j 


= E 




X ^(-*i)l{xi,>o,E.«,,,<*^, i^^._^^^,^=o}_ 




/•oo 

= dt/(t)//(e-''^*l|x,>o})E 
Jo 


L ^ iex,ti>-t '' ^ 






X hm E 
i<:-s.oo 


exp — A 




0}i 


)J 



where the first equality is based on the values of A and the second one holds since Poisson point 
measures over disjoint sets are independent. We will calculate the terms respectively. 
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First we have 

^(e"'^^* 1|X,>0}) = /«(1{X.>0} - (1 - e-"^'))- 
Using the expression of /i, we have 

/■oo 

^(1{X,>0}) = /3Qo(^t > 0) + / n{dx) P,(Xt > 0) 

Jo 

POO 

= /3c{t)+ / n{dx){l-P^{Xt = 0)) 
Jo 



Pc{t) + / n{dx) (1 - e-^-^W) = F(c(t)), 

JO 



and 



/•CX3 

/i(l - e-^^0 = /3Qo(l - e-^^0 + / ^(^a;) P.(l - e"''^*) 

JO 
/•oo 

= I3vt{r]) + / n(dx) (1 - e-^'^*^'')) = F{vt{ri)). 
Jo 

Second, using Lemma [3?T] we get 



E 



exp( -7 Y^ XU^ 

iei,U>-t 



exp 



F{vs{-f))ds). 



Finally we see that 



Um E 

K^oo 



exp 



exp 



(-A Yl iXU-u+Kl^xu^o})) 



dsl 



{s>o}/"(l-e~^''=l 



{Xs+t=0}j 



exp r - y ds l|,>o}/u(l - e-^^= PxA^t = 0))") 

exp(-|dsl|,>o}M(l-e-(^+^W)^0^ 
(isF(t;,(A + c(t))) y 



= exp 

where we use the exponential formula for the Poisson point measure in the first equality and 
the Markov property of X in the second one. 

Putting all the calculations together we obtain the result. □ 

It is then straightforward to derive the distribution of the TMRCA A. 
Corollary 4.1. The distribution of A is given by 

F{A < t) = exp ( - / F{c{s)) dsj = E[e-'^(*)^o], 

and A has density Ja with respect to the Lebesgue measure given by 

(/■oo N 

- J F{c{s)) ds 
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Proof. Using Theorem 14.11 and (j2.6p , we see that the first equality holds easily. The second one 
is immediate. □ 

We see that in this general case the expression of the distribution of A is invariant compared 
with [8] . The next result is also a direct consequence of Theorem 14.11 

Corollary 4.2. Conditionally on A, the random variables Z , Z , Z are independent. 

We also derive from Theorem 14.11 the distribution and the mean of the population size just 
before MRCA. As can be seen from below that the expression for the Laplace transform is the 
same as that of [8]. 

Corollary 4.3. Let t > 0. Then 

(4.1) E[.-^-^lA = t]=^i^^ ani E[Z-|A = .] = £#». 

^ ' L I J ]E[g-c(t)Zo] t I J '0(c(i)) 

Proof. This is a direct consequence of Theorems 14.11 and Corollary 14.11 □ 

We can further deduce that conditionally on {A = t}, the distribution of Z converges to the 
distribution of Zq as t — > oo. 

Another application of Theorem 14.11 we call the bottleneck effect, is that the size of the 
population just before MRCA is stochastically smaller than that of the current population. 
Note that this inequality does not hold in the almost surely sense in general. The proof is the 
same as that of [H]. 

Corollary 4.4. For all z > and t > 0, we have F{Z^ < z\A = t) > F(Zo < z). Hence the 
population size Z is stochastically smaller than Zq, that is F[Z < z) > F{Zq < z), for all 
z >0. In particular we have E,[Z \A] < K[Zo]. 

Remark 4.1. Instead of considering the size of the population just before MRCA, we consider 
the size at MRCA, Z^, which is given as Z^ = Z^ + ^iei ^o^{ti=-A}- We don't take into 
account the contribution of i € J72 since for those we have Xq = 0. Similar calculations as those 
of Theorem 14.11 show that for A,t > 0, 



E[e-'^^+\A = t]=E[e-^^ \A = t]- 



F(c(t)) 

If F (0) = oo, then linit^-oo E[e~ + |j4 = t] = 0, which means that conditionally on {A = t}, Z^ 
is likely to be very large, as i — )■ cxj. We can interpret it as this: a clan is born at time —t and 
it survives up to time 0, if t is large enough, it is likely to have a large initial size. Therefore, 
Z^ is not stochastically smaller than Zq in general. 

5. The number of oldest families 

In this section we will consider the number of families in the oldest clan alive at time 0. It 
is equivalent as that of individuals involved in the last coalescent event of the genealogical tree. 
We will use the family representation in this section. 

The number of oldest families alive at time (excluding the immortal individual) is defined 
as: 



(5-1) ^^ = Ev=-*„xi„>o} = Ei{ 



^{A=-tj, xi^ >o} ~ Z^ ^{A=-tj, Qy-tj}- 



Obviously A^"^ > 1. In particular when /3 > and the measure n = 0, we have A^"^ = 1. 
The following theorem gives the joint distribution of A,N and Zq. 
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Theorem 5.1. Let < a < 1. For any non-negative measurable function f , we have 



E[a^^e-^^«/(^)] 



Proof, For i G X, set 



dsf{s)exp( - / F{vr{X))dr- / F{cir))dr 

\ Jo Js 

X (f{c{s)) - F((l - a)c{s) + aVsiX)) 



J* 



Ji^i, if i e /, 
{i}, ifieJ2- 



For / non-negative measurable, we have 

E[a^^e-^^«/(A)] 



E 



g-^Efcei,tfe 




dsf{s)E 


Q~^ EfeGI,tfc<0 ^-tfe l{tfe>-s} 




(/3Qo[ae-^ 


""^ l{x.>o}] + / n(dx)E, 

JO 





where the first equality is based on the decomposition of A, the second on splitting tj. into three 
parts: t^ > s^t^ < s and t^ = s, and XligJa ^x^i^X) is a Poisson point measure with intensity 
xQo^dX) under P. We will calculate the terms separately. 

By Lemma |3.H we have 



E 



exp( -A Yl ^-t,Mt,>-s} 
fcex,tfc<o 



exp 



and 



El 

fcGX 



{tk<-s,X'l^>0} 



I = exp 



F{vr{X))dr 



F{cir))dr]. 



The next equation is obtained by splitting the terms into two parts: 



Q^ieJa ^{xi>o} 



Xi>0} 



E, 



^iGJs {Xl>0} 



expr-A^Xf")] -pJ^Xf=0 



The first part is calculated as follows: 



Ex 



.S.6./3 1{ 



X|>0} g-^EjsJg^s 



exp f - xQo [1 - al|x,>o} e ^^^ ] j 

exp { - xQo[X, > 0] + xaQo [l{x,>o} e"^^^ ] j 

exp ( - X [(1 - a)c{s) + aUs(A)] j . 
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The second part is 

-xQoCl-e-^^s) 






e 

/ A — vnn 

ieJa 

= lim e-^"^(^) = e-^'^^") 

A— >-oo 



For the last part, we see that 

Qo[e-^^^ l{x.>o}] = Qo[^s > 0] - Qo[l - e"^^^] = c{s) - u,(A). 
Putting all the calculations together we see that 

/■CO / pS /"CO > 

K[a^^ e-^^'' f{A)] = / dsf{s)exp{ - / F{vr{X))dr- / F{c{r))dr 

Jo \ Jo Js J 

X [f{c{s)) - F((l - a)c{s) + av,{\))\ . 
This finishes the proof. □ 

Using the density of A, the following corollary is immediate. 

Corollary 5.1. For < a < 1, A,t > 0, we have 

^l N^ -xZoiA ,1 F{c{t))-F{{l-a)c{t) + avt{X)) f f^ 

E[a e o|^ = ij= F(c(t)) exp I - / F{vr{X))dr 

and 

.... ^r7V-|, ,1 F{cit))-F{{l-a)c{t)) F((l - a)c(t)) 



The next corollary is direct from Corollary 15.1 
Corollary 5.2. VKe have for n > 1, 



.A _,,_._, ,,.+ic(t)"FW(c(t)) 



F[N'' = n\A = t] = {-ly 



n! F(c(i)) ■ 



Then E[iV^|^ = t] = ^j.°t|f|*^ G [0, oo]. /n addition if F' (0+) < oo, i/ie function t ^ 
^[N \A = t] is non-increasing. 

Notice that if we let F{t) = ct", c > and < a < 1, then the conditional distribution will 
not depend on the CB-process but only on the immigration structure. 

6. The number of ancestors at a fixed time 

In this section we will consider the number of ancestors Mg at time —s of the current pop- 
ulation living at time and how fast it tends to infinity. To answer this question we need to 
introduce the genealogy of the families which is a richer structure studied in [9l [lOl [21 [1] . 



6.1. Genealogy of CB-process. The construction developed by Duquesne and Le Gall [9l[T0] 
for (sub)-critical CB-process is well known. Results in [lOj is restated in the framework of the 
measured rooted real trees, see [2]. We will follow Section 2 in [1]. 
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6.1.1. Real tree framework. A metric space {T,d) is a real tree if the following two properties 
hold for every s,t £ T, 

• (unique geodesic) There is a unique isometric map fs^t from [0, d{s, t)] into T such that 
fs,tiO) = s and /s,t(d(s,t)) = t. 

• (no loop) If q is a continuous injective map from [0,1] into T such that g(0) = s and 
q{l) = t, we have q{[0, 1]) = /,,i([0, d{s, t)]). 

A rooted real tree is a real tree (T, d) with a distinguished vertex called the root. Denote such 
a tree by {T,d,^). If s,t G T, we will note [[s,t| the range of the isometric map fs^t described 
above. We also denote [s,i[[= ls,t} \ {t}. 

If x € T, the degree of x, n(x), is the number of connected components of the set T\ {x}. 
The set of leaves is defined as Lf (T) = {x £ T\{0}, n{x) = 1}. The skeleton of T is the set of 
points in the tree that are not leaves: Sk(T) = T\Lf(T). 

For every x E T, |0, xj is interpreted as the ancestral line of vertex x in the tree. If x, y G T, 
there exists a unique z £ T, called the Most Recent Common Ancestor (MRCA) of x and y, 
such that [0, x] n [0, y| = |0, zj. Then the root can be seen as the ancestor of all the population 
in the tree. We shall call the height of x, h{x), the distance d(0,x) to the root. The function 
X !-)> h{x) is continuous on T, and we define the height of T by Hmax{T) = sup^,g7-/i(x). 



6.1.2. Measured rooted real trees. We will denote by T the set of the measured rooted real trees 
(T, d, 0, m) where (T, d, 0) is a locally compact rooted real tree and m is a locally finite measure 
on T. We may simply write T in case of no confusion. 

Let T G T. For a > 0, we set T{a) = {x € T, d(0, x) = a} for the level set at height a, and 
T^a{T) = {x € T, d(0,x) < a] for the truncated tree T up to level a. We consider iTaiT) with 
the root 0, d'^'^^'' and m^"^'' are the restrictions of d and m to TTaiT). Let {T^'°, k € IC) he the 
connected components of T\ TTa{T). Denote by 0^ the MRCA of all the vertices of 7"'^'°. Set 
T^ = T^'° u {0fc} which is a real tree rooted at point 0^ with mass measure m' defined as the 
restriction of m' to T^ . We will consider the point measure on T x T: 

k&K. 



6.1.3. Excursion measure of Levy tree. Recall that ^ is a (sub)-critical branching mechanism. 
There exists a c-finite measure (or an excursion measure of Levy tree) N[d7^ on T, with the 
following properties: 

(i) (Height). Va > 0, N[F^ax(r) > a] = c{a). 

(ii) (Mass measure). The mass measure m' is supported on Lf(T), N[(i7^-a.e. 
(iii) (Local time). There exists a T-measure valued process {i°',a > 0) cadlag for the weak 
topology on finite measure on T such that N[(i7^-a.e.: 



nJ{dx) = / t{dx) da, 
Jo 



e^ = 0, inf{a > 0;r = 0} = sup{a > 0;r / 0} = i?max(T) and for every fixed a > 0, 
N[d7l-a.e.: 

• The measure £"" is supported on T{a) . 
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• We have for every bounded continuous function (p on T: 



lim ^ j <^(^)l{H^ax(r')>e}A/'a^.(c^a;, dT'). 



Under N, the process {{(."•, 1), a > 0) is distributed as X under Qq. 
(iv) (Branching property). For every a > 0, the conditional distribution of the point measure 
J\fJ{dx, dT') under N[(iT|-ffmax(7') > a], given vra(T), is that of a Poisson point measure 
on T{a) X T with intensity r((ix)N[dr']. 

In order to simphfty notations, we will identify X with {{i"", 1), a > 0) as well as Qo with N. 
We give a definition for the number of ancestors. 

Definition 6.1. The number of ancestors at time a of the population living at time b is the 
number of subtrees above level a which reach level b > a: 



Ra,b{T) - 2^ l{//„,„,(r'=)>6-a}- 

fce/c 
6.2. Genealogy of stationary CBI-process. We use l\3.2\i to construct the genealogy of Zq. 

• Conditionally on A^O) let Ni{dt,dT) = J2jeJi \t,T^)idt,dT) be a Poisson point measure 
with intensity v{dt)Q,Q{dT) with v{dt) = Ylii^i'''i^ti{dt). 

• Let N2{dt,dT) = "^j^j^^t ,T3)idt,dT) be a Poisson point measure independent of 
No,Ni with intensity (3dtQo{dT). 

We will write X^ for i"'{T-^) for j € J^. Thus notation (13. 2p is still consistent with the previous 
sections, "^j^j ^{t-^T^) allows to code the genealogy of the family of Zq in the Levy tree sense. 

Let s > 0, we will consider the number of ancestors at time — s of the current population 
living at time 0, that is, 

3&J 



6.3. Asymptotic for the number of ancestors. First we present the following theorem. 
Theorem 6.1. The conditional joint distribution of Mg and Zq is: for rj, X > 0,s > 0, 

In particular, conditionally on Z^g, Mg is distributed as a Poisson random variable with 
parameter c{s)Z-s. 
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Proof. For any r], X > 0, we have 



E[e^ 



-bZ.^-rjMs-XZo] 



E 



exp 



exp 



F{vr{X))dr 



{-s<tj<0}' 






xE 



exp 



E i{*.<-}(^ ^-.-t, + ^^~-*„-*. (r^) + Axi, 



jej 



exp 



/ F(t;^(A)) dr\ exp ( - / daF[Qo(l - Q-^^--^^-,--rs{r)-\x^^, ^j \ 



where in the first equality we use the fact that the Poisson point measures over disjoint sets are 
independent, in the second one we use Lemma l3. II and we use an immediately generalization of 
Lemma l3.1l to genealogies in the third equality. Using branching property we have 



Qc 



1 



-hXa-riRa,a + siJ)~\XaA 



Qo 
Qo 



1 _ e-^<i{6+Qo[l-cxp(-»?l{^,,^^^(7-)>3}-AXs)]} 



Since X^ = on {Rmax{T) < s}, we have l-e-''^{«— (^)>''>-^^'' = {l-e-'')l{^>,} + e-'^{l- 
Q-XXs\ Then we deduce that 



Qo 



1 



-bXa~riRa,a + a(T)-XXa + s 



Qo[i 



-XXa 



Vai>^ 



with A = 6 + (1 — e ^)c{s) + e ^ iJsW- Then we get the result. 



D 



Intuitively Mg counts the number of excursions of the height process at time —s above level 
s. Similar result as that of Duquesne and Le Gall |1U] can be deduced here. 



Corollary 6.1. The following convergence holds: 

M, 



lim 



Zq, a.s. 



s^o c{s) 

7. The MRCA age process and the zero set of the CBLprocess 

In this section we will deal with the MRCA age process {At , t G M) and the zero set of the 
CBI-process Z = {t £ E., Zt = 0} by using the sample path decomposition of the CBI-process. 

For the clan decomposition of the CBI-process shown in Section 3, we have used the Pois- 
son point process Ns{dt,dX) = '^i^x\u,x^){dt,dX) with intensity dtiJ,{dX), where fi{dX) = 
/3Qo(d-^) + /(-Q QQ-, n{dx)Fx{dX). The corresponding Poisson point process which is in charge of 
the birth time and duration time is N4{dt,dC) = Y2iel ^{ti,C')i'^^^ "^0 with intensity dt^{C, € •), 
where /x(C > t) = F{c{t)). 

7.1. The MRCA age process. Define the left leaning wedge with apex at (i,r) by 

A(i, r) := {{u, v),u < t andn + v > r}, 

which is the set of points that give birth before t and is still alive at time r. We can thus define 
the MRCA age process {At,t G M) as follows: 

At:=t- inf{s : 3C > 0, such that (s, C) G {{U, Q) : i e 1} H A(t, 0)}. 
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The strong Markov property of the Poisson point processes A'4 imphes that {At,t G M) is a time 
homogeneous Markov process. The transition probabihties of the MRCA age process can be 
proved along the same hnes as that of part (a) of Theorem 1.1 in Evans and Ralph |11] . 

Theorem 7.1. The transition probabilities of the Markov process {At,t G M) are seperated into 
two parts: 

• for < y < X -\-t, the continuous part 

nAs,. e dy\A^ = X) = fl - ^^^1^) e-/r-W"))-F(c(,))d,, 

• otherwise, a single atom ¥{As+t = x + t\As = x) = p(5^)) ■ 

7.2. The zero set of the CBI-process. From the clan decomposition, we have 

• either Zt 7^ 0, for t £ [ti,ti + Ci),i € /; 

• or Zi / 0, for t G (tj, ti + Ci),i e J2. 

Then we have 

{t,Zt = 0} = R\ \J[t^,ti + Ci)U \J{ti,t, + Q) 

Then we can derive that the zero set of the CBI-process is a random renewal set with the 
following equation: 

Z = R\\J{ti,ti + Q). 

iGl 

To see why this is always true, we only need to focus on those points that belong to M \ 
UiGx(*«'*« + Ci) but not to {t,Zt = 0}. These points actually are the left accumulation points 
of {t,Zt = 0}. Indeed suppose that this is not true. Let {ti} be such point. Then for any ij, 
there exists e^ > such that [ti — ei,ti) n {t,Zt = 0} = 0. This is impossible since after taking 
the closure we obtain that [ti — ei,ti]r\ Z = $; while ti G Z. 

We derive the following theorem. 

Theorem 7.2. (1) Z = % if and only if f^ exp {f^ F{c{u))du) dt = 00. 

(2) If Jq exp ( Jj F{c{u))du) dt < 00, the radom set Z has a positive Lebesgue measure a.s. 
if and only if Jq°° -jjrw dt < 00; the random set Z is the union of the closed intervals of 
positive lengths if and only if f3 = and n{dx) < 00. 

Proof. We have known that /i(C > t) = F{c{t)). We can derive directly from Corollary 5 in [15] 
that the first assertion holds. 



With a slightly modification of Proposition 1 in [15] or Propositon 1.22 in [M], we can obtain 
that Z has a positive Lesbesgue measure if and only if 

/•OO /'OO /"OO 

/ tfi{C£dt)= fj,{C>t)dt= F{c{t))dt < 00. 

Jo Jo Jo 

Letting r = c{t) in the above equation yields J^ ^44 dt < 00. In order to derive the last 
statement, we use Corollary 2 in |15j . which requires //((" G dt) to be finite, i.e. F{c{t)) is finite, 
as t ^- 0+. Since we have c{t) -^ 00 as t ^- 0+, we need /3 = and n((0, cxd)) < 00. □ 

Easy calculation gives a simple example that ipiu) = 2l36u + /3u'^ and F{u) = 2(3u satisfying 
condition (1), and Z = $. Another example is given in [13] in stable case. 
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